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ABSTRACT 


The NOAA-7, GOES-5 and GOES-6 VISSR/VAS solar channels have been calibrated 
for the periods from October 1983 through January 1985 (NOAA-7, GOES-6) and from 
October 1983 through July 1984 (GOES-5). The White Sands Monument area in New 
Mexico, whose relectance properties are well known, and space are used as calibration 
targets. The shortwave reflected terrestrial irradiance that is measured at satellite altitude 
is computed using a fairly accurate radiative transfer model which accounts for multiple 
scattering and bidirectional eflects (Tanrd et al., 1979). The groimd target reflectance and 
relevant characteristics of the overlying atmosphere are estimated from climatological data 
(ozone amoimt, aerosol size-frequency distribution and refractive index) and observations 
at the nearest meterorological sites (water vapor amount, visibility). The approach is 
believed to produce accmacies of 8 to 13% depending on the channel considered. The 
results obtained for the 15-month period do not indicate any drift in the solar channels of 
both types of instrument. The gain changes, about 15% of the mean values, are largely 
attributed to inhomogeneities of the ground tzurget (shading eflects due to the presence of 
dunes). No systematic eflect of the normalization procedure applied by NOAA to the raw 
VISSR/VAS data is detected. There is some evidence that the GOES-5 solar channels 
gradually deteriorated from March 1984 until the failure of the satellite in July 1984. 
Comparisons of gains determined in orbit with those before launch show that the NOAA-7 
solar channels read higher by about 15%. The disparities cannot be explained by model 
errors and must have occurred before the time period analyzed here. 


1. Introduction 


The Visible Infrared Spin Scan Radiometer/Vertical Atmospheric Sounder (VISSR/VAS) 
and the Advanced Very High Resolution Radiometer (AVHRR) aboard the meteorologi- 
cal satellites operated by the National Oceanic and Atmospheric Administration (NO A A) 
were primarily intended to provide frequent, high resolution, visible and infrared images 
for monitoring weather events and temperature stratifications over large areas. Sucessfully 
applied to cloud detection, the data in the visible and near-infrared bands have also proved 
suitable for other environmental applications, in particular surface albedo and vegetation 
index mappings {e.g., Rockwood and Cox, 1978; Tarpley ct al., 1984; Tucker et al., 1984), 
and radiation budget studies [t.g., Stephens et al., 1981; Gautier, 1982, 1986). For these 
applications, the need for accurate calibration of the involved satellite data is easily per- 
ceived, and can hardly be overemphasized. To illustrate this. Fig. 1 shows the impact of 
calibration errors on selected geophysical quantities currently derived from satellite data, 
namely the normalized vegetation index, the surface albedo and the downweu'd solax irra- 
diance at the surface. The impact is particularly dramatic on the normalized vegetation 
index (Fig. la). For a green biomass level of 660 gm~*, a 20% error in the calibration of 
the AVHRR visible (channel 1, 0.58-0.68 ^m) and near-infrared (channel 2, 0.73-1.10 jum) 
bands translates into a typical 61% error in the normalized vegetation index. The error 
rapidly increases with decreasing biomass level, and reaches 114% for a 260 gm~^ biomass 
level. In the case of the surface albedo (Fig. lb), a 20% error in the calibration of the 
VISSR/VAS visible band (0.52-0.73 /xm) introduces a similar error in the surface albedo, 
typically ranging from 19 to 34% between the extreme cases of white gypsum sand and 
clear sea water. Finally, it is clear from Fig. Ic that for highly reflecting clouds (t.c., albedo 
>0.5) even small calibration errors in the VISSR/VAS visible band yield important solar 
irradiance errors. When monitoring such quantities over long-term periods, for instance 
in climate studies, the smallest fluctuations are often determinant and must, therefore, be 
resolved. If the satellite data used in the calculations sure not accurately calibrated, not 
only will small changes not be detected, but the results may sdso turn out unreliable smd 
misleading [t.g., comparisons of different data sets become nonsensical). 
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Fig. 1. Sensitivity of vegetation index, surface insolation and ground 
albedo to calibration errors. 





Careful preflight laboratory calibrations of the VISSR/VAS and AVHRR visible chan- 
nels were performed, respectively, by the Santa Barbara Research Center (Staff Members, 
1980) and by ITT Aerospace (Lauritson tt al., 1979; Koczor, personal communication). 
Transfer functions were established to calculate satellite level irradiances from the teleme- 
tered digital sensor signab. These functions, however, may not be valid for the data ac- 
quired in operations because the orbital satellite environment differs from the steady state 
laboratory environment and inflight changes occur, as observed with other radiometers. A 
serious problem is the possibility of dark current from the detectors and zero offsets in the 
signal conditioning circuitry. Degradation of the optics with time is also not uncommon. 
It is therefore important to maintain a check-of-calibration for the visible channels while 
the instruments operate in orbit. 

Unfortunately, the AVHRR contains no calibration capabilities other than space view- 
ing for chatnnels 1 amd 2. Onboaurd facilities exist for the VISSR/VAS visible channel (t.e., 
am opticad subaissembly with reduced size aperture for directly viewing the sun), but they 
are not being utilized because of imcertainties in the chauracteristics of the optical elements. 

Several substitutive methods of calibration can yet be envisioned. One can utilize 
independent amd coincident measurements from a suitably equipped high-flying aircraft 
{t.g., Smith amd VonderHamr, 1980; Kriebel, 1981; Hovis ct ai., 1985) or by a calibrated 
satellite-bome instrument with similar spectral chairau:teristics (e.g., Smith amd Loramger, 
1977; Brooks et al., 1984). Another possibility, perhaps the simplest way, as pointed out by 
Coulson amd Jacobowitz (1972), is to use aureas on the eaurth’s surface as calibration targets 
{c>g., Koepke, 1982). This last method is adopted in the present study amd is applied to 
NOAA-7 AVHRR, GOES-5 and GOES-6 VISSR/VAS shortwavelength chamnels. In this 
paper, we report on the results obtadned for the period from October 1983 to JamuauT* 1985. 


2 



2. Method of Calibration 

The first step is to select a suitable calibration target on the earth’s surface. Once the 
target is selected, the intensity of the radiation that is directed to space from the target 
to the satellite b calculated from the optically acting parameters of the target and of the 
intervening atmosphere, and is related to the sensor digital output. The procedure is then 
repeated periodically to monitor eventual time changes in the behavior of the instruments. 

The GOES-5 and GOES-6 VISSR/VAS vbible channeb are actually comprbed of eight 
independent detectors, similar in spectral response (Figs. 2a and 2b) and sharing the same 
telescope and scan assembly. The analog signal from the detectors is digitized into 6-bit 
counts before its transmission to the grotmd, but two bits are added to each count at the 
NO A A Command and Data Administration (CD A) station in Wallops bland (t.e., users 
have access to 8-bit coded data with, however, only 6-bit accuracy). While the detectors 
respond linearly to incident irradiance, the digitbing scale b not linear, but purposely 
adjusted to keep the signed-to-nobe ratio linear with respect to digital coimt, t.e., the 
digital count is proportional to the square root of the volt response (Fig. 3). Thus, the 
irradiance sensed by each VISSR/VAS vbible detector varies as 

E = GX^ + 1 (1) 

where E is the irradiance of the target at count value X, and G and I are calibration 
constants. It b important to emphasbe that thb relationship, currently adopted by in- 
vestigators, b only approximate. At count value 48 on the 6-bit scale, for instance, a 7% 
error in b introduced when using the best fit of Fig. 3 for the analog-digital transfer 
function. 

As opposed to GOES-5 and GOES-6 VISSR/VAS visible channels, for which the spec- 
tral response of the detectors defines the optical bandwidth, the wavelength range of 
NO AA-7 AVHRR channels 1 and 2 (Fig. 2c) is controlled by a filter. The voltage output 
of the detectors (one per channel), converted into 10-bit counts aboard the spacecraft, is 
linear as a function of sensed irradizince (Lauritson et al., 1979), i.t., 

E=GX + 1 (2) 

describes the relationship between counts and irradiances. 

To determine the “constants” G and J, two calibration points (two JC, E pairs) are 
needed. The White Semds Monument area in New Mexico (hereafter referred to as White 
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Fig. 2. Spectral response of the solar channels. 






Fig. 3. Analog-digital VISSR/VAS transfer function. 




Sands), located at 32®52'N, 106^17'W and 1.22 km altitude, which is highly reflective, 
provides a suitable high calibration point. First, the area is Izirge, approximately 30 km 
in diameter, providing multiple sampling points and minimizing edge and environmental 
effects. Second, its reflectance properties are well-known, stable in time and fairly well 
approximate those of a Lambert siirface, reducing bidirectional effects and simplifying the 
radiative transfer modeling. Finally, favorable atmospheric conditions mostly prevail in 
that region (minimum cloudiness, relatively low aerosol content). The area, however, is 
not uniform horizontally. A high water table keeps the very low paxts in a permanently 
wet state, modifying its reflectance. Pzirt of the area is covered by dunes, and their shading 
is difficult to compensate for, thereby restricting the observations to conditions of small 
solar and viewing zenith angles. The selection of small zenith angles is also important, as 
we shall see, to reduce uncertainties resulting from tmknown atmospheric characteristics. 
The second calibration point, at practically zero intensity value, is provided by space. 

As evidenced in Figs. 2a and 2b, the spectral response of GOES-5 and GOES-6 
VISSR/VAS visible detectors may differ in the longwavelength half power point by as 
much as 0.1 /xm, which suggests that the calibration be effected separately for each de- 
tector. However, an average calibration of the eight detectors is accomplished for each 
satellite, using average count data and average spectral responses (also shown in Figs. 2a 
and 2b). This procedure is justifled because the data which are made available to users 
by NOAA CDA have been processed to eliminate displeasing stripes that arise from dif- 
ferences among the detectors. The output of the eight detectors is normalized to one of 
the detectors, and the normalization is readjusted every couple of weeks. Moreover, in 
some cases, different detectors are used as the normalizing detector. In other words, the 
data calibrated one week may be out of calibration the following week. Consequently, 
our vicarioTis calibration of the GOES-5 and GOES-6 VISSR/VAS visible channels will 
be performed every two weeks and, for convenience, we shall proceed in the same way for 
NOAA-7 AVHRR channels 1 and 2. 

3. Radiative Transfer Model 
a. Model 

In the spectral region of interest, 0.5-1. 1 /xm, the sensed energy originates almost 
entirely from solar radiation which propagates through the atmosphere, is reflected by the 
surface target and is Anally transmitted to the space-borne sensor. The thermal radiation 
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emitted by the atmosphere itself is negligible in comparison with the solzir component, so 
it will be ignored. 

Solar radiation in the 0.5-1. 1 fua wavelength range is scattered by air molecules and 
aerosols, and is absorbed primarily by ozone, water vapor, oxygen and aerosols. Scatter- 
ing and absorption processes complicatedly interact, but fortunately, gaseous absorption 
can be treated separately (Deschamps et al., 1983; Tanre et al., 1985). Qualitatively, 
the arguments are the following. First, ozone is located at altitudes where molecules 
and aerosols are already very rarified, i.t., the ozone layer is traversed almost without 
scattering. Ozone, therefore, simply acts at reducing the radiation signal by the direct 
transmittance (t.e., along a direct path] corresponding to the double passage across the 
absorbing layer. Second, water vapor absorption occurs at wavelengths greater than 0.85 
^m, where molecular scattering is insignificant and where only aerosol scattering subsists. 
Since (l) aerosol scattering is mainly primary and secondary at those wavelengths, and (2) 
the aerosol phase function presents a pronounced forward peak, the photons practically 
follow a direct path through the atmosphere. To account for water vapor absorption, it 
is therefore sufficient to afiect the signal that would be observed if the atmosphere were 
not absorbing by the transmittance along the direct path from the sun to the target and 
from the target to the satellite. Part of the solar radiation incident at the top of the atmo- 
sphere is in fact scattered back into space before any surface refiection takes place. The 
separation of scattering and absorption processes is more difficult to justify in this case, 
but remains a reasonable approximation since absorption by atmospheric gases is generally 
weak or very localized in the 0.5-1. 1 tim range. This also warrants the same treatment for 
oxygen absorption, which occurs in a very narrow band near 0.76 /im, although molecular 
scattering cannot be considered vmimportant at this wavelength. 

Thus, the irradiance measured by the space-borne sensor viewing the surface target 
can be written with a good degree of approximation as 

a f FxEQxcosdQtgx{Ooi6) pxiM,6o,0,<l>)dX (3) 

Jo 

where A is wavelength, 0q and 6 are respectively the solar and viewing zenith angles, (j> 
is the relative azimuth angle, Fx is the spectral response of the detectors or filters, ^oA 
is the extraterrestrial solar irradiance (in Wm“*/xm“^), tgx is the transmittance due to 
absorbing gases, and p\is an equivalent reflectance of the target M. 
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The radiation geometry is depicted in Fig. 4. The equivalent reflectance p\ is deflned 


as 


a*(M e 6 !h\ — 
p^{M,eoA<t>) - Eoxcose 


( 4 ) 


where is the irradiance that would be sensed if the air molecules were non-absorbing. In 

other words, p\ takes into account aerosol and molecular scattering and aerosol absorption, 
but excludes gaseous absorption. The study deals with monochromatic quantities, but 
subscript A will be omitted afterward to simplify the notation. 

The equivalent reflectance p* is modeled according to Tanrd et al., (1979): 


p*{6oJ,4>) ^Pa{M,eoA<l>) + i—(p)i^M:0oA<f>)r 

[Ap{M,0o,0,4>) + Bp{M,0oA<t>) + C{p){M,0oA*t>)] 


( 5 ) 


where Pa is the intrinsic atmospheric reflectance (t.c., the contribution of the signal back- 
scattered into space without surface reflection), r is the spherical albedo of the atmosphere, 
p is the bidirectional reflectance of the target, and "p and (p) are average angular and envi- 
ronmental reflectances. In Eq. 5, A, B and C are development coefficients characterizing 
the relative importance (or weight) of reflectances p, p and (p). Detailed expressions for 
A, B and C can be foimd in Tanrd et d., (1979). 

The reflectance p corresponds to photons going diffusely from the sun to the surface, 
and then from the surface directly into space, t.e., 


C"" f 0, ^ - <f>) EJ {00, 0\ <!>•) cos 0> sm0> d0>d<l>> 

0 ^ 0 


p{M,0o,0,<i>) — 


f f ^ Ej{0Q,0',<IP)cos0’sm0'd0'd<IP 

JQ Jq 


( 6 ) 


where EJ is the downward diffuse irradiance at the surface for p = 0. 

The reflectance (p) represents photons that have sustained reflections from the area 
adjacent to the target, and appears in Eq. 5 because the reflectance of this surrounding 
background a priori differs from that of the target. It is deflned as 


/ +00 f +00 

I COB 0oEqp(M\ 0', 0, 4) t{0, x', yQdx'dy' 
^pnm,oo,fr,9) = 


f f ^ Ej {00, 0',<tP) COB 0>Bxa.0’d^'d<fP 

JQ Jo 


( 7 ) 


where t{0,z^,y') denotes the efficiency of a point M*{^,j/) with horizontal coordinates x* 
and y'. 
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Fig. 4. Radiation geometry. 


b. Solar irradianct 


The solar irradiance incident at the top of the atmosphere is taken from Neckel and 
Labs (1984) and corrected for earth-sun distance variations during the course of the year. 
The total irradiance in the different channels, as obtained by convolving Neckel and Labs’ 
(1984) data with the respective response functions, is given in Table 1. Clearly, more energy 
flux is available in the spectral dom^ of the VISSR/VAS channels, for the reason that 
their peak of transmission is closer to the wavelength of maximum solar emission, and also 
because their bandwidth is relatively large. For AVHRR channel 2, the large bzmdwidth 
of the Alter partly compensates for the smaller energy flux emitted by the sun in the 
near-infrrured. In order to facilitate the conversion of total irradiances into monochromatic 
irradiances, we also listed in Table 1 the values of the integral over wavelength of the 
different filter functions. 

c. Atmospheric functions 

To evaluate the gaseous transmittance tg, water vapor, oxygen and ozone are considered 
separately. Following Tanr^ ct al., (1985), the individual transmittances are parameterized 
as a function of air mass and vertically integrated amount of absorber. The statistical band 
models of Goody (1964) for water vapor and Malkmus (1983) for ozone and oxygen are used 
in the calculations. The spectral resolution is 20 cm~^, which allows a precise description 
of the phenomena (t.e., absorption strength and overlapping of bands). To account for 
the dependence of the absorption coefficients on temperature and pressure, the Curtis- 
Godson approximation is employed with the US62 standzurd atmosphere of McClatchey 
(1971). When absorption by several gases takes place, tg is computed as the product of 
each separate transmittance. 

The atmospheric functions pg, r and the coefficients A, B and C (which also depend 
on the characteristics of the atmosphere) were calculated by Tanre ct al., (1985) for dif- 
ferent experimental conditions and various atmospheric models. Fairly accurate analytical 
formulas were established, which give />a, r. A, B and C as a function of viewing geome- 
try and optical thickness of the atmosphere (or turbidity). As an alternative, the surface 
visibility, a parameter routinely measured at all weather stations, may be used instead 
of the turbidity. This approach is appealing because turbidity is only measured at a few 
sites, and is selected for the present study. One should note, however, that the correspon- 
dence between visibility and turbidity is only a rough inverse proportionality (thick layers 
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of aerosols may exist aloft with a very clear atmosphere at the stirface). Because of the 
lack of precise knowledge of the aerosol scattering properties at White Sands (i.e., size- 
frequency distribution, refractive index), we opted, among diverse aerosol models available, 
for the continental model of the International Radiation Commission (WCP, 1983). Al- 
though this model may not represent what actually occurs at White Sands, it certainly 
remains useful in setting limits to what might be reasonably expected. 

d. Surface reflectance 

The bidirectional reflectance of White Sands p was measured in laboratory by Wal- 
raven and Coulson (1972) for dry and wet conditions. In-ettu measxirements were also 
carried out by Hovis (1966), but they are not applicable to the present study. Hovis 
(1966) employed a spectrophotometer that measured monochromatic energy in a direction 
normal to the stirface and, consequently, his results do not represent the range of smgular 
conditions encountered in the real world (the White Sands reflectance is a relatively strong 
function of the angle of incidence). In fact, the energy collected by the integrating sphere 
of the instrument b reflected in all directions into the hemisphere in front of the sample 
surface, i.e., angular effects are lost. 

Fig. 4 summarizes the main features reported by Walraven and Coulson (1972). For 
dry samd, the spectral reflectance curve (Fig. 5a) is characterized by a rapid increase to 
0.6 pm, followed by a slight dip around 0.7 pm, continued by a barely perceptible increase 
beyond 0.8 pm. For wet sand, the entire reflectance curve is displaced toward lower values 
by about 0.1, or 15%. The slope is steeper below 0.6 ^m, and the depression at 0.7 
pm is more pronoimced. Above 0.8 pm, as opposed to dry sand, the curve exhibits a 
negative declivity caused by water vapor absorption in the 0.94 pm band. The shape of 
the curves below 0.6 pm is consistent with the general finding that, for sandy textured 
soils, the reflectance sharply drops in magnitude at wavelengths smaller than the average 
particle size. In interpreting Fig. 4a, however, one should bear in mind that the portion 
of the curves extending beyond 0.9 pm results from a linear extrapolation of Walraven 
and Coulson’s (1972) data. This questionable procedure is at least consistent with Hovis’ 
(1966) measurements, which exhibit a fairly linear comportment in the 0.8-1.0 pm domain. 

Important angulzir effects at 0.7 are clearly seen in Figs. 4b and 4c. Variations 
with the solar zenith angle (or equivalently with the viewing zenith angle because of the 
reciprocity principle) are particularly significant for dry sand (Fig. 5b), i.e., 15% when 
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Fig. 5. Reflectance of White Sands zis a function of wavelength, 
solar and viewing zenith angles, and relative azimuth angle 
(after Walraven and Coulson, 1972). 





passing from 0 to 60° for the specified conditions. The refiectance curve for wet sajid is 
characterized by a nearly constant slope, and the overall effect is more moderate (only a 5% 
change in the 0-60° range). The behavior is similar at other wavelengths (not shown here): 
generally, an isotropic distribution is observed at small zenith angles, all the more so as the 
wavelength decreases. Almost no dependence on the relative azimuth angle is discerned for 
dry sand (Fig. 5c). This is in agreement with the previous statement, because the zenith 
angles considered in Fig. 4c are small. At larger zenith angles, a specular component 
superimposes the diffuse contribution, its importance increasing as the principal plane is 
approached (i.e., ^ = 180°). Even for solar and viewing zenith angles of 15° and 30°, 
however, the wet sand data (Fig. 5c) shows definite evidence of a speculax component. 

Inspection of Fig. 5 clearly reveals the strong infiuence of surface moisture upon White 
Sands refiectance, which therefore ought to be carefully taJcen into account. Unfortunately, 
soil moisture is not currently measured at White Sands and, anyhow, the usefulness of 
punctual or isolated soil moisture data, even recorded when the satellite views the target, 
is debatable. A correction is therefore estimated with the help of gravimetric measure- 
ments performed by Williamson (1975) during a one-year period. These measurements 
are reproduced in Fig. 6 for two sites, one located in the fiat portion of the area and the 
other in the white gypsum sand dimes. Clearly distinguishable in Fig. 5 is a seasonal de- 
pendence, with minimum values during summer (about 10%) and maximum values during 
winter (typically 30% in the fiats and 20% in the dunes). Short-term variations are also 
evident, and can be of importance {e.g., 20% over a few day period around Julian day 200 
in the fiats). The yearly averages are 21.4 ± 5.3(l<r)% in the fiats and 17.8 ± 3.2(1<t)% 
in the dunes. Since the standard deviations for the entire one-year period translate into 
surface refiectance changes of less than 0.6%, we concluded that the yeaxly averages would 
adequately represent the actual moisture conditions at White Sands. We axe aware, how- 
ever, that this procedure is somewhat arbitrary, because of the high spatial and temporal 
variability of surface moisture. Also, the measurements of Williamson (1975), even though 
they correspond to clear or partly cloudy skies without rain, were taken only near noon 
(local), and we do not know about the atmospheric events preceding these meaLSurements'. 
A better way to proceed would certainly consist of directly measuring the refiectance of 
White Sands at the very moment of the satellite pass, but such a practice is beyond 
the scope of this study. Although the distinction between fiats and dunes in not always 
apparent in the satellite data, we decided, for consistency, to restrict our analysis to the 
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course of the year (after Williamson, 1975). 



dunes area. In the simulations, consequently, the surface moisture is fixed at the yezurly 
average value of 17.8%. Since Walraven and Coulson’s (1974) data indicate that, in some 
degree, the refiectance of white gypsum sand veiries linearly when the sand dries from a 
saturated state, the moisture effect at the intermediate value of 17.8% is simply computed 
by linear interpolation. 

e. Average reflectances 

The average reflectances p and {p) are laborious to evaluate, because their calculation 
{via Eqs. 6 and 7) involves angular integration of the diffuse irradiance EJ and requires 
precise knowledge of the spatial variations of the target reflectance. Fortunately, the 
difficulty can be obviated. In order to do this, let us first examine the relative importance 
of development coefficients A, B and C (Fig. 7). One notices that, in the wavelength 
range of interest, A Isurgely surpasses B and C in magnitude for clear atmospheres (t.e., 
for visibilities greater than 23 km) and not-too-horizontal sightings. In other words, the 
contribution of the signal directly reflected by the ground, which naturally involves the 
true reflectance of the teirget p, is highly preponderant for sufficiently clear atmospheres. 
Second, the spherical albedo of the atmosphere r does not exceed the value of 0.2, even 
for considerable aerosol loadings (Fig. 8). Third, p and (p) do not notably depart from 
p. As shown in Fig. 9, (p — p) /p lies below 5% in most situations. It decreases with 
increasing solsur zenith angle and, for the conditions specified in the figure, diminishes to 
zero at solax zenith angles of 40-50°. Beyond these values, the difference is negative (t.e., 
p < p) . This feature simply reflects the influence of the specular component of the White 
Sands reflecteince, which is important at large solar zenith angles, even in the direction 
perpendicular to the principal plane. In the computation of p, EJ was expressed as the 
sum of an isotropic component and a Dirac function (see Tanre et al., 1983), which makes 
p directly proportional to p. For given zenith angles, (p — p)/p decreases with increasing 
relative azimuth angles, but moderatly in the case of Fig. 9 because the dependence of p 
on is small at ^o=15° and ^=30° (relatively small zenith angles). The variations with 
the viewing zenith angle, almost negligible at \—0.5p,m in Fig. 9, are notable at A=0.9^m, 
and become more important 2 is ^ increase. Because EJ appears in both the 

numerator and denominator of (6), the dependence of p — and a fortiori of {p — p)fp — upon 
the atmospheric visibility is partly smoothed out. 
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Fig. 9. Departure from angular average reflectance p of p for differ- 
ent wavelengths, model atmospheres and geometrical con- 
ditions. 
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The angular variations of {{p} —p)/{p) (Fig. 10) are almost nonexistent, a direct conse- 
quence of the large horizontal «ctent of the White Sands area. The atmospheric visibility, 
a parameter of minor importance at A = 0.9 ^m, significantly affects (/?) and, hence, 
{(p) — p)l{p) at A = 0.5 pm, because scattering processes are more operative at smaller 
wavelengths. Because the reflectance of the White Sands environment is lower than that of 
White Sands, (p) is smaller than p, and the magnitude of the difference increases with the 
visibility (with fewer aerosols, more photons originating from the environment are scat- 
tered in the direction of observation). This b clearly seen in Fig. 10: ((p) — p)/ (p) decreases 
from about —0.5% to —5% when the visibility is increased from 5 to 23 km. To evaluate 
(p). White Sands was assumed to be a circular target, 30 km in diameter, surrounded 
by an homogeneous background of reflectance p^ = 0.3 and, following Deschamps et al. 
(1983), (p) was written as a linear combination of p and pj. 

The preceding considerations, and the fact that (p — p)/p and ((p) — p)/(p) are of 
opposite signs, strongly incline one to replace p and (p) in Eq. 5 by p. In so doing, the 
modeling of the apparent reflectance p* b drastically simplified and reduces to the case 
of an homogeneous Lambertian surface. The validity of the procedure is substantiated 
in Fig. 11, which shows the resulting error on p*. Clearly, the error mainly depends on 
atmospheric visibility and zenith angles, and does not exceed ±1% in most conditions. 
The apparent reflectance p* b either underestimated or overestimated, depending on the 
visibility and on the radiation geometry. This interestingly indicates that, by neglecting 
the influence of the environment, the selection of clear atmospheres might not be optimum, 
particularly for those channels whose peak of transmission is located near 0.5/xm. In the 
presense of hazy atmospheres, however, the errors introduced by our insufficient knowledge 
of aerosol optical properties are increased. Fig. 11 abo indicates that Ap*/p* is generally 
reduced at small zenith angles. Such angles axe also favored to minimize errors arising from 
uncertainties in the characteristics of the atmosphere overlying the target (the atmospheric 
path is smaller). In many cases, however, and thb b especially true for a sun-synchronous 
satellite, it will not be possible to meet the conditions (visibility, radiation geometry) that 
correspond to the most desirable compromise. 
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Fig. 11. Error introduced in the modeling of p* by replacing p and 
{p) with p. Positive values indicate underestimation. 
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Fig. 12. Error in the modeling of p* in the case of Lambertian and 
homogeneous ground. The comparisons are made with ex- 
act calculations involving the successive orders of scattering 
method. Positive values indicate underestimation. 









/. Accuracy 

In order to assess the accuracy of our modeling of the irradiance E measured by the 
satellite-borne sensors, errors in the parameterization of g 2 iseous absorption and on the 
model of Tanre et al., (1979) should be considered in addition of the errors that we just 
discussed above. According to Tanre et al., (1986), the error on tg is smaller than 1%, 
except at grazing incidence and observations angles {0q,$ > 80®). The model of Tanre et 
al. (1979) has not been validated against in situ measurements, even in the simple case 
of an homogeneous and Lambertian ground (a difficult task anyhow). Comparisons with 
exact calculations involving the successive order of scattering method indicate, however, 
that for turbid atmospheres with visibilities of about 23 km, the model underestimates p* 
by less than 1% in the wavelength range of interest (Fig. 12). At A=0.5/iim , this error 
happens to partly compensate for the error introduced by our approximation of p and (p), 
but not at A=0.9/im (see above). In fact, by restricting our observations to reasonably 
clear conditions (visibility of ~23 km), the absolute value of the overall error expected on 
our modeling of E should not exceed 2%. 

4. Data 

a. Satellite data sets 

The VISSR/VAS data used for the calibration were provided by the Space Science and 
Engineering Center of the University of Wisconsin, Madison in the form of full resolution 
(«.€., ~1.5 km at the latitude of White Sands), 6-bit precision {i.e., 1/64) 256x256 pixel 
images coded on an 8-bit scale. The images were acquired at approximately 15 day inter- 
vals, 2 ind cover the periods from October 1983 through January 1985 for GOES-6 eind from 
October 1983 through July 1984 for GOES-5 {i.e., until the failure of the satellite). For 
each selected day, the data set consists of three consecutive (i.e., half hourly separated) 
clear sky images of White Sands and one image of space, all taken near the local noon. 

The NOAA-7 AVHRR data were chosen among orbit passes received in real time and 
archived at the Satellite Oceemography Facility of the Scripps Institution of Oceanography, 
San Diego. They also cover the period from October 1983 through January 1985. White 
Sands generally falls within the acquisition range of the ground station, but not always, 
and is often located near the east edge of the 2900 km-wide swath. The peissage over the 
antenna site occurs at about 16:00 local time, which is relatively late in the afternoon and. 
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Fig. 13. Examples of White Sands images in the AVHRR solar 
channels. 
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Fig. 14. Same as Fig. 13, but VISSR/VAS instead of AVHRR 






consequently, not optimum for the calibration, especially during winter (too low sun above 
the horizon). The direct readout High Resolution Picture Transmission data telemetered 
to the station are 10-bit in precision 1/1024) and the ground resolution at nadir is 
1.1 km. 

Figs. 13 and 14 give examples of White Sands images in the VISSR/VAS and AVHRR 
shortwavelength channels, respectively. White Sands is easily distinguishable from the 
surrounding environment by its high brightness (white tones). The shape of the area 
appears somewhat distorted from one image to the other, a direct consequence of the 
different viewing geometries. The demarcation between the flats and the more reflective 
dunes, not discernible on the VISSR/VAS images (Fig. 13) is clearly apparent on the 
AVHRR images (Fig. 14). The contrast, in fact, depends on the geometrical conditions 
of observation. The VISSR/VAS images exhibit stripes, which illustrate the differences 
among the detectors (see section 2). This suggests that the destrip ing procedure applied 
by NO A A CD A might not be completely efficient, since the corresponding space images 
(not shown here) were relatively free of stripes. 

Sub-images of the White Sands dunes were extracted from the original satellite data 
sets. Since no geometric corrections were applied to the data, the number of pixels con- 
tained in the sub-images is not constant, but depends on the viewing geometry (we seek 
to include as many pixels as possible to improve the statistics of the calibration results). 
Since our objective for the VISSR/VAS is to perform an average calibration of the eight 
shortwavelength detectors, eight consecutive lines form the sub-images. The space counts 
at the beginning of the AVHRR scan lines that correspond to the White Sands viewing 
and the entire 256x256 pixel VISSR/VAS images of space constitute the zero intensity 
data set used in the computations. 

Average White Seinds and space counts and their respective standard deviations are 
presented in Figs. 15 and 16 for AVHRR and VISSR/VAS shortwavelength channels, re- 
spectively. Over the study period, the average space count is practically constant for 
AVHRR channels 1 and 2, and does not significantly change for the VISSR/VAS short- 
wavelength channels. In particular, no appreciable shift or trend can be detected. The 
average White Sands count, on the contrary, strongly varies with time. The low frequency 
flucuations reflect those of the solar radiance at the top of the atmosphere since, for each 
satellite, the observations were taken at approximately the same time during the entire 
study period. In some instances, abrupt changes are observed at about two week intervals 
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Fig. 15. White Sands and space AVHRR brightness counts for the 
study period. 
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Fig, 16. Same as Fig. 15, but VISSR/VAS instead of AVHRR. 
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{c.g., between February 2 and Februzury 8, 1984 for GOES-5) or from one image to the 
next during the same day {e.g., on December 30, 1984 for GOES-6). The spatial (i.e., 
horizontal) standard deviation is generally small (2 to 6 counts) for AVHRR channels 1 
and 2, but often exceeds 10 counts for the VISSR/VAS shortwavelength channels, and may 
result from surface moisture variations, inhomogeneities in the surface composition, local 
changes in the characteristics of the atmosphere, intrinsic instrumental noise, differences 
among the detectors (for GOES-5 and GOES-6) and instrument changes. 

6. Viewing geometry 

The angles characterizing the incidence and observation directions are displayed in 
Figs. 17, 18 and 19 for GOES-5, GOES-6, and NOAA-7, respectively. 

During the study period, the solar zenith angle varied from 10 to 60® for GOES-5 and 
GOES-6, and from 45 to 80® for NOAA-7. The low frequency fluctuations exhibit a sinu- 
soidal form, which is intimately related to the annual cycle of the sun. The NOAA-7 solar 
zenith angles are generally large, which makes their selection for an accurate calibration 
questionable. But we are limited to the afternoon NOAA-7 passes, and they occur when 
the sun is already low. 

The viewing zenith angles from the geostationary altitude remain fairly constant (».e., 
~50® for GOES-5 and ~49® for GOES-6) omtil GOES-5 failed in July 1984. GOES-6, 
stationed over the Equatorial Paciflc at 135® W at the moment of the failure, was then 
maneuvered over the United States to the longitude of 85®W, and the viewing zenith angle 
changed to about 38®. The NOAA-7 viewing zenith angle is subjected to more important, 
somewhat erratic changes (Fig. 17), which in fact depend on the disteince of the satellite 
subtrack to White Sands (the closer the satellite passes the area, the smaller is the viewing 
angle). 

Successful capture by the tracking antenna generally implies that NOAA-7 peiss west- 
ward of the target. Since, on the other hand, the data are acquired when the sun is in the 
western quadrant, the sun and the satellite are viewed by the target from about the same 
azimuthal direction, thereby leading to small relative azimuth amgles. On one occasion 
(t.e., November 30, 1984), however, the satellite passed close to, but east of White Sands 
and a relative azimuth angle of nearly 158® was observed. For GOES-5 and GOES-6, the 
relative azimuth angle rapidly varies from one hour to the next (Figs. 18 and 19). This is 
simply due to the fact that the observations are taken near the local noon, t.c., near the 
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Fig. 17. NOAA-7 viewing geometry and atmospheric characteristics 
at the moment of satellite overpass. 
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Fig. 18. Same as Fig. 17, but GOES-5 instead of NOAA-7. 
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Fig. 19. Same els Fig. 17, but GOES-6 instead of NOAA-7. 




maximum solar elevation during the day. The amplitude of the fluctuations is larger during 
su mm er, because the sun is higher at noon during that season. Note also in Figs. 18 and 
19 that before the failure of GOES-5, the GOES-5 eind GOES-6 relative azimuth angles 
vary in opposite direction during the day. The reeison for this is that the two satellites 
were viewing White Sands from directions located respectively eeist and west of the target. 
Indeed, after the failure of GOES-5, the same type of variations as those observed for 
GOES-5 axe encountered for GOES-6, except that the magnitude of the relative azimuth 
angle is smaller, a consequence of the larger GOES-6 azimuth angle (the GOES-6 subpoint 
on the earth is closer to the White Sands meridian). 

c. Meteorological parameters 

The vertically integrated ozone amounts were estimated from the meteorological atlas 
of London et al., (1976). At the latitude of White Sands, the values range from 0.27 cm.atm 
in autumn to 0.32 cm.atm in spring. The vertically integrated water vapor amounts were 
obtained from radiosonde observations at El Paso (31®48'N, 106®24'W), a meteorological 
station located at nearly the same altitude and longitude as White Sands. No corrections 
for the latitudinal distance between El Paso and White Sands (about 120 km) were elfected. 
This procedure is roughly correct because the atmospheric perturbances usually travel 
west in that region. The atmospheric visibility is that visually observed at meteorological 
station HMN, a station very close (a few km) to the calibration site. Unfortunately, a 
standard value of 25 km was reported continuously during the entire study period. The 
actual visibility might in fact be much larger in many cases; values of 100 km or more have 
been frequently observed (Slater, personal communication). Our deflcient knowledge of the 
visibility at White Sands is indeed a potential source of error for the calibration, which can 
be estimated (see below) but not accounted for. At White Sands, the atmospheric pressure 
is considerably reduced because of the surface altitude. This was carefully accounted for in 
the computations of the Rayleigh optical thickness. Atmospheric pressure also affects the 
absorption strength of the intervening gases as in the case of temperature, but the error 
that is introduced when using the US62 standard proflles of pressure and temperature is 
neglected because of its secondary importance. 

5. Results 

The satellite irradiances computed for White Sands, with the data and the model 
discussed above, exe presented in Fig. 20 for NOAA-7 and in Fig. 21 for GOES-5 and 
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Fig. 20. Modeled White Sands AVHRR irradiance for the study 
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Fig. 21. Same as Fig. 20, but VISSR/VAS instead of AVHRR. 




GOES-6. They correspond to the brightness counts of Figs. 15 and 16, respectively. Large 
variations on the seasonal time scale are clearly present in Figs. 20 and 21, of the order 
of 80 and 60Wm"^ for AVHRR channels 1 and 2, respectively, and of 80Wm~* for the 
VISSR/VAS visible channels. For consecutive times during a given day, the VISSR/VAS 
irradiances vary by only a few Wm“* because the sun zenith angle does not change much 
during the one-and-one-half hour period of the observations. 

The modeled irradiances fairly well reproduce the seasonal or even monthly variations 
of the counts, but not those existing at smaller time scales, especially from one hour to 
the next in the case of GOES-5 and GOES-6. This is partly due to our method of specify- 
ing the atmospheric parameters, which does not account for changes occurring at hourly 
time scales. Also, the short-term variations of the counts may result from instrument 
modifications, which indeed are not considered in the method. 

To evaluate the accuracy of the simulated irradiances, uncertainties of 30% on the 
anisotropy factor of the aerosol phase function, 30% on the atmospheric visibility, 10% on 
the ozone amount, 20% on the water vapor amount, 50% on the soil humidity and 1% on the 
surface refiectance (due to other processes than soil humidity) were assumed. Uncertainties 
on other aerosol parameters, such as their vertical distribution, are not significant to the 
first order. The 10% uncertainty on the ozone amount corresponds to the daily fluctuations 
around the geographical monthly mean (London et al., 1976). The 20% uncertainty on the 
water vapor amount is typically that of radiosonde observations. The 50% uncertainty on 
the soil humidity is representative of the annual variations reported by Williamson (1975). 
The 1% uncertainty on the surface reflectance accounts for errors in the measurements 
of Walraven and Coulson (1972). The 30% uncertainty in the asymmetry parameter is 
cheiracteristic of differences between the extreme cases of maritime and urban aerosols. 
Finally, the 30% uncertainty on the atmospheric visibility is somewhat arbitrary, because 
the visibility at the calibration site is poorly known. But the major problem with the 
visibility undoubtedly arises from the standard value of 25 km consistently reported by 
meteorological station HMN. Is this value representative of the average conditions at the 
calibration site? If not, and recent observations tend to indicate that the visibility at 
White Sands is often larger (>100 km), the modeled irradiances are consequently biased 
(».c., underestimated) by several percents, more precisely by 2.5% if the average visibility 
is 100 km. Thus, except for the visibility, the uncertainties estimated for the atmospheric 
and surface parameters are reasonable, and would most likely maximize the overall error 
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Table 2. Uncertainties on model paraoneters and 
corresponding errors in simulated irradiances. 


A£7/JE; 

(%) 

AP/P 

P {%) VISSR/VAS AVHRR/1 AVHRR/2 


Anisotropy Factor 

30 

0.5 

0.7 

0.6 

Visibility 

30 

0.9 

1.1 

1.0 

Water Vapor Amount 

20 

0.2 

0.2 

1.2 

Ozone Amount 

10 

0.6 

0.7 

<0.1 


Surface Reflectance 
Soil Humidity 


1 

50 


1.5 

2.2 


1.4 

1.9 


1.3 

2.1 


on the modeled irradieinces. On top of these uncertainties, one should not forget the 
uncertainty on the radiative tranfer model itself (about 2% for most situations encountered 
in the study; see Section 3e). When introduced together in the model, the uncertainties 
discussed above yield the errors shown in Figs. 20 Eind 21, computed for all individual 
cases. These errors range from 8 to 13%, and might explain the inability of the method to 
retrieve the short-term variability of the brightness counts. 

To apprehend what are the most critical parameters for the calibration, the different 
uncertainties were considered separately around average conditions for the atmosphere. 
The results are presented in Table 2. Clearly, the most contributing errors are those on 
the soil humidity, the surface reflectance zmd on the atmospheric visibility. The influence 
of the anisotropy factor is relatively limited, because the effects of this parameter on pa on 
the one hand eind on A, B and C on the other hand (see Eq. 5) partly compensate. For 
the ozone amount, our treatment based on climatological data appears to be sufficient. 
For the water vapor amount, if the utilization of radiosonde observations is satisfactory 
for the VISSR/VAS visible channels and for AVHRR channel 1, it yields an error of about 
1.2% for AVHRR channel 2. This error, as most of the other errors, can be reduced, for 
instance by using optical instruments on site (c.g., pyrheliometer for the turbidity, Voltz 
photometer for the water vapor amount). This would require, however, increased means 
and effort, and will be discussed in the final section as part of a global strategy to control 
the quality of the data acquired in orbit. 

From the White Sands and space average counts (Figs. 15 and 16), the modeled irradi- 
ances (Figs. 20 eind 21), and assuming that the space irradiance is null, average calibration 
coefficients were calculated. This was done by comparing the modeled irradiances with 
the corresponding AVHRR 10-bit counts (Fig. 22) or VISSR/VAS 8-bit counts squared 
(Fig. 23). The correlation coefficient is better than 0.995 at the 99% confidence level for 
both types of instrument. Assuming that the calibration points should fit a straight line of 
the form (1) or (2), the best estimates of the slope and intercept were computed by means 
of the principle of maximum likelihood, and are given in Table 3. In the computations, 
however, the linear regressions were forced toward the space points, because the irradiance 
of space (equal to zero) is not subjected to uncertainty and, further, the time variability 
of the space counts is small. If this constraint were not introduced, the results would have 
changed significantly for GOES-5 and GOES-6 because the VISSR/VAS calibration points 
are not well distributed, yielding less accurate calibration coefficients. The residual errors 
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. 22. Plots of modeled White Sands AVHRR irradiance versus 
brightness counts. Statistics of best linear fits are specified 
on the figure. 
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Fig. 23. Same as Fig. 22, but VISSR/VAS instead of AVHRR. 





Table S. Average post-launch calibration coefficients. 


S atellite/Sensor 

Gain 

(Wm“^/coTintt 
or Wm“^/count^t) 

Intercept 

GOES-5 VISSR/VAS 

5.080x10-3 ±0.069x10-3 

-2.5 ±1.8 

GOES-6 VISSR/VAS 

5.127x10-3 ±0.064x10-3 

-3.2 ±1.5 

NOAA-7 AVHRR/1 

0.2137 ±0.0020 

-7.5 ±0.4 

NOAA-7 AVHRR/2 

0.3154 ±0.0041 

-11.5 ±0.7 


fAVHRR 

iviSSR/VAS 


axe small, i.e., about 5% of the mean irradiance value for the VISSR/VAS visible channels 
(Fig. 22), and only 3% and 4% respectively for AVHRR channels 1 and 2 (Fig. 23). This 
provides a good idea of the level of stability of the calibration coefficients, and tends to 
indicate that our assessment of the errors on the modeled irradiances is pessimistic, but 
the residual errors do not include eventual biases. 

6. Discussion 

One of the goals of the calibration is to assess whether the instruments are experiencing 
changes {e.g., drift) as a result of performance degradation in flight. This can be verifled 
by exeimining the evolution with time of the calibration coefficients, and by comparing 
the results with those obtained in laboratory before launch. In the following, the cases of 
AVHRR channels 1 and 2 and of the VISSR/VAS visible channels are discussed separately. 

o. AVHRR channels 1 and 2 

Fig. 24 depicts the evolution with time of the AVHRR channels 1 and 2 calibration 
coefficients. In the flgure, the error bars associated to individual points represent the 
differences due to spatial variations in the White Sands counts. Clearly, there is no evidence 
of a systematic drift in any of the computed coefficients and no dramatic change occured 
during the 15-month period of the observations. The gains fluctuate somewhat erratically 
throughout the study period, by 11 and 15% (lo) of the mean values for channels 1 and 2, 
respectively. The larger excursion of channel 2 gain likely results from the higher sensitivity 
of this channel to water vapor (the effect of uncertainties on the water vapor amount is 
incre^lsed). These fluctuations, however, do not appear to be signiflcant at the level of 
accuracy of our calibration; the errors in the modeled irradiances may easily account for 
the observed variability. 

It is interesting to note that the gains are highly correlated to the viewing zenith 
angle. The correlation coefficients are —0.77 (channel 1) and —0.83 (channel 2), signiflcant 
at the 99% confldence level. This intriguing feature cannot be attributed to errors in 
the parameterization of gaseous absorption, or even to the fact that absorption-scattering 
interactions were neglected, because absorption processes act quite differently in the two 
channels, which would not produce similar correlation coefficients. Also, when examining 
the model errors and their dependence on the viewing geometry (Figs. 9 and 10), there is 
no evidence that an increase in viewing zenith angle would be reflected in a decrease of the 
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Fig. 24. Time series of AVHRR calibration coefficients determined 
in orbit. Pre-launch values are represented by dashed lines. 



gains. On the contrary, one would conclude that p* is slighty more overestimated at higher 
viewing zenith angles, yielding a positive correlation. One possible explanation is linked to 
surface reflectance inhomogeneities. In fact, shading effects are enhanced when the sun is 
low {recall that 6 > 50® for most conditions), which reduces the apparent target reflectance 
at small viewing zenith angles because incidence and observation directions have similar 
azimuth angles. In other words, p* is more overestimated as 0 decreases because more 
shaded areas occupy the instrument field of view; these areas are implicitly assumed to be 
illuminated by the stm in the model. The effect is almost impossible to quantify, however, 
because this would require precise knowledge of the shape and geometry of the dunes and 
other spatial inhomogeneities; but it might realistically explain the irregular variations of 
the gains observed in Fig. 24. 

When compared to prelaunch gains, the average gains obtained in the present study 
read higher by about 14% for channel 1 and 15% for channel 2 (Fig. 24). The differences 
are significant, and cannot be explained by model errors: these errors do not contain biases 
of such magnitudes. In fact, uncertainties in the atmospheric visibility likely introduce a 
negative bias (see section 5), which would aggravate the disparities. Hovis (1982) pointed 
out several problems in the technique applied by ITT Aerospace to calibrate the AVHRR 
solar channels, from the standpoint of the calibration in the laboratory with a sphere and 
from the fact that there still is no agreement on the value of the solar irradiance in the 
bands of interest. We used Neckel and Labs’ (1984) solar data instead of Thekaekara’s 
(1974) tables, but this cannot quantitatively explain the discrepancies. A likely possibility 
is therefore that a degradation of the solar channels happened after launch. Since Fig. 24 
indicates the occurrence of little or no drift, the solar channels might have deteriorated 
suddenly. It would certainly be worthwhile to check whether a degratation took place as 
a result of launch constraints. 

So far we have only analyzed changes in the channels 1 and 2 gains. Comparison of 
intercepts determined in orbit with those before launch (Fig. 24) show that the post-launch 
values are lower by 2 Wm”^ (channel 1) and 3Wm“* (channel 2). The time variations 
reflect those of the gains, but naturally are inversely correlated since two-point calibrations 
were performed. We recall here that the zero offsets (in counts) were found practically 
constant for the two channels, which is evidence of the stability of the system. 
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b. VISSR/VAS visible channels 

In order to assess the stability of the VISSR/VAS solax channels with time, average 
calibration coefficients were computed for each selected day of the study period. For a 
given day, two-point calibrations were performed with the observation of space and each 
of the three White Sands observations, and the results were averaged to provide the daily 
value. Fig. 25 shows the evolution with time of the gains and intercepts obtained in that 
way. The error bars in the figure represent the variability of the coefficients during the 
one-and-one-half period of the observations. The gain cheinges, about 12% (GOES-5) and 
15% (GOES-6) of the meam values, are not erratic as those of AVHRR channels 1 and 
2. Still, they remain within the uncertainty domain of our calibration (c.g., model errors, 
transfer function uncertainties and the limitations due to the 6-bit digitization). We note, 
however, that the VISSR/VAS viewing zenith angles, unlike the AVHRR ones, did not 
vary much throughout the study period, except for GOES-6 immediately after the failure 
of GOES-5 (t.e., when GOES-6 wets maneuvered to its new position). The solar zenith 
angle, on the contrary, varied noticeably and a correlation coefficient of 0.66, significant 
at the 99% confidence level, is computed between the GOES-6 gain and the solar zenith 
angle. In fact, at fixed or minimally varying viewing zenith angles, p* and, consequently, 
the gains axe more overestimated as increases because, here again, the instruments view 
more shaded axesis. This might be the cause of at least part of the variability in the GOES- 
6 gain. Such a high correlation, however, is not observed for GOES-5. On the contrary, 
from March 1984 until the failure of the satellite, period during which 6o decreases, an 
increeise in gain is noted (Fig. 25). This reversed trend might therefore be associated to 
changes in the system that eventually lead to its collapse. 

No systematic effect of the normalization procedure can be detected. However, an 
abrupt increase in the GOES-5 gain (i.e., a loss of sensitivity) of about 17% of the mean 
value is observed between February 2 and Februeiry 18, 1984, which cannot be explained 
by shading effects or inhomogeneities in the surface reflectance. The fact that on March 3 
the gain recovers its value of February 18 might be the manifestation of the normalization 
procedure or other intervention at the ground station. Note that the phenomenon occurred 
just before the solar channels began to lose sensitivity, which might therefore indicate that 
something drastic and irreversible happened to the system in orbit at that time. Important 
variations axe also observed during some days, for instance on October 4, 1983 for GOES-5 
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Fig. 25. Same as Fig. 24, but VISSR/VAS instead of AVHRR. 



ajid December 30, 1984 for GOES-6, and correspond to those cases for which the White 
Sands brightness counts considerably vary over a half-hour interval. Because the changes do 
not last, they cannot originate from the normalization procedure. A possible explanation 
is an error at the ground station that is immediately corrected for. 

From the pre-launch calibration curves provided by the Santa Barbara Research Center 
and the actual scanner temperatures communicated by the NO A A GOES support branch, 
the pre-launch gains and intercepts were computed and are displayed in Fig. 25 against 
the post-launch values. The sc 2 mner temperatures consequently vary from about 31° in 
winter to 10° in summer, which translates in a significant increase in the sensitivity of 
the GOES-5 detectors (8%) and a slight decrease (<1%) in that of the GOES-6 detectors. 
Deviations of up to 4% may exist, however, from these average figures between certain 
detectors. The seasonal changes in the heating effect are not retrieved at all in the post- 
launch values of Fig. 25; they are probably masked by the more important variations in the 
apparent refiectance discussed above. The thermal effect on the intercepts is obviously not 
visible. These, in other respects, depart from the pre-launch values by — lWm“* (GOES- 
5) and — 3Wm“* (GOES-6), which generally represent a negligible contribution to the 
signal measured by the instruments. In a word. Our analysis thus indicates a fairly good 
agreement between pre- and post-laimch calibration coefficients, the discrepancies being 
largely attributed to spatial inhomogeneities in the surface refiectance, except for GOES-5 
from March 1984 until the failure of the satellite period during which the solar channels 
plausibly deteriorated. 

7. Siumnary and Recommendations 

The potential of measurements from space of the solar radiation reflected from the 
eeurth-atmosphere system for environmental applications has been extensively illustrated 
in the recent literature. Vegetation index, ground and cloud albedoes, surface insolation 
and radiation budget at the top of the atmosphere represent a non-exhaustive list of 
parameters that are discussed in terms of satellite observations in the visible portion of 
the solax spectrum. Most of the satellite instruments to be used in determining these 
parameters, unfortimately, axe incompletely chaxacterized and calibrated, except those on 
the Earth Budget Radiation Experiment. Such is the case of the VISSR/VAS and AVHRR 
radiometers onboard the meteorological satellites operated by NOAA, which do not have 
active inflight calibration capabilities for the solar channels. 
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A simple, probably lezist expensive, way to routinely calibrate the solax channels of 
these instruments while they operate in orbit is to use areais on the earth surface as cali- 
bration targets. This technique was applied in the present study to the AVHRR onboard 
NOAA-7 Eind to the VISSR/VAS onboard GOES-5 and GOES-6. The White Sands Mon- 
ument area in New Mexico, of well-known bidirectional reflectance, was used as reference 
standard. Observations taken when the sensors were pointing the deep space served to fix 
the offset points on the count- versus-irradiance curves. The satellite level irradiances were 
computed from the ground target characteristics and the optically acting parameters of 
the overlying atmosphere using a fairly accurate radiative transfer model, which accounts 
for multiple scattering and bidirectional effects. Atmospheric properties were determined 
either from climatological data (ozone amount, aerosol phase function) or from observa- 
tions at the nearest meteorological site (water vapor amount, visibility). The calibration 
accuracy is believed to be on the order of 8 to 13%, depending on the viewing geometry, 
surface and atmospheric conditions and on the channel considered. 

The results obtained for the 15-month period analyzed indicate the occurence of no 
significant drift in the solar channels of both types of instrument. The gain changes 
(about 15% of the mean values) are largely attributed to spatial inhomogeneities of the 
calibration site (shading effects due to the presence of dunes) that are not accounted for in 
the modeling. There is some evidence, however, that the GOES-5 solar channels gradually 
deteriorated from March 1984 until the failure of the satellite in July 1984. Comparisons of 
gains determined in orbit with those before launch show that the AVHRR solar channels 
read higher by about 15% on the average. This apparent loss of sensitivity cannot be 
explained by model errors, whose bias component tends to aggravate the discrepancies. 

The accuracy of the calibration could certainly be improved by measuring m situ and 
at the moment of the satellite overpass all the necessary parameters. This would require, 
however, increased means and effort. In other respects, the non-negligible uncertainties 
introduced by the inhomogeneity of the target might be avoided by performing the cali- 
bration in the flat portion of the White Sands eirea. A detailed validation of the radiative 
transfer model is also necessary. In order to do this, a complete data set on atmospheric 
properties and stirface reflectance (including radiation measurements) has been recently 
acquired for the Yuma desert in Arizona. The Yuma desert, in fact, appezirs to be a 
near-homogeneous, very suitable site for the calibration. Its use with White Sands would 
allow a better verification of the linearity of the detectors with respect to incident energy, 
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and more confidence would be placed in the results. In ajiy case, the calibration of the 
GOES VISSR/VAS and NOAA AVHRR solar channels should be adressed in the context 
of a global and concerted strategy {e.g., see lAMAP, 1981; WCP, 1986). In view of the 
requirements for many applications (in particular climate studies) , the only way to ensure 
that the observed variability in the satellite-derived parameters is the result of changes 
on eeirth and not in the instruments is to calibrate the data acquired in operations on a 
continuous basis. As demonstrated in the present study, this can be achieved with our 
simple and relatively inexpensive method. Nevertheless, episodic (».c., one or twice during 
the life of individual satellites) but intensive process-oriented campaigns that involve the 
utilization of high-fiying aircraft should be conducted, even if they are onerous, for they 
constitute the optimal way of assessing the data quality. 
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